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We give a full analysis of the auto- and cross-correlations between the Stokes parameters of 
the cosmic microwave background. In particular, we derive the windowing function for an antenna 
with Gaussian response in polarization experiment, and construct correlation function estimators 
corrected for instrumental noise. They are applied to calculate the signal to noise ratios for future 
anisotropy and polarization measurements. While the small-angular-scale anisotropy-polarization 
correlation would be likely detected by the MAP satellite, the detection of electric and magnetic po- 
larization would require higher experimental sensitivity. For large-angular-scale measurements such 
as the being planned SPOrt/ISS, the expected signal to noise ratio for polarization is greater than 
. . . one only for reionized models with high reionization redshifts, and the ratio is less for anisotropy- 

polarization correlation. Correlation and covariance matrices for likelihood analyses of ground-based 
and satellite data are also given. 

vn 
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The detection of the large-angle anisotropy of the cosmic microwave background (CMB) by the COBE DMR 
experiment Q provided important evidence of large-scale spacetime inhomogenities. Since then, a dozen of small- 
, scale anisotropy measurements have hinted that the Doppler peak resulting from acoustic oscillations of the baryon- 
I ' photon plasma on the last scattering surface seems to be present CMB measurements gain an advantage over 
, other traditional observations due to the fact that the small CMB fluctuations can be well treated as linear, while 
' the low-redshift universe is in a non-linear regime. It is now well established that CMB temperature anisotropics 
T-H ' are genuine imprint of the early universe, which could potentially be used to determine to a high precision virtually 
C . all cosmological parameters of interest. It has been estimated that a number of cosmological parameters can be 
\ determined with standard errors of 10% or better by the upcoming NASA MAP satellite Furthermore, the future 
f~| Planck Surveyor CMB mission would have capability of observing the early universe about 100 times better than 

map. 

Q , At this point, we should explore as much information as possible besides the temperature anisotropy contained in 
' the relic photons. Anisotropic radiation possessing a non-zero quadrupole moment acquires a net linear polarization 
' when it is scattered with electrons via Thomson scattering (also see Eq. (6) of Ref. (|]). When the photons begin to 
d decouple from the matter on the last scattering surface and develope a quadrupole anisotropy via Sachs- Wolfe effect [|, 
^ linear polarization is created from scatterings with free electrons near the last scattering surface. Studies have shown 
• i that on small angular scales the rms polarization, in a standard universe, is a few percents of the rms anisotropy, while 
/\ '■ the large-scale polarization is insignificant |7j . In models with early reionization, the large-scale polarization is greatly 
enhanced, to a few percents level, but the small-scale anisotropy is suppressed significantly ||,^. Therefore, CMB 
polarization would provide a valuable complementary information to the anisotropy measurements. In addition, the 
anisotropy-polarization cross correlation offers a test of physics on the last scattering surface, as well as a possibility 
of distinguishing the scalar and tensor perturbations However, all of these polarization calculations have relied 

on a small-angle approximation, which may not be valid when a large sky-coverage is considered. 

As such, full-sky analyses of the polarization have been performed ]ri|~p^ . It was found that there are modifications 
to low multipole moments [l < 30) of the polarization power spectra, where the tensor contribution dominates over 
the scalar contributio n pd^ jll] . More importantly, rotationally invariant power spectra of the Stokes parameters have 
been constructed P4-|15||. In particular, one of them is a parity-odd magnetic polarization spectrum, which vanishes 
for scalar-induced polarization, thereby allowing one to make a model-independent identification of non-scalar (i.e. 
vector or tensor) perturbations (also see Ref. |l6| ). Recently it was shown that magnetic polarization would be a 
strong discriminator between defect and inflation models pT^JlS]. Also, a new physically transparent formalism based 
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on the total angular momentum representation was proposed which simplifies the radiative transport 

problem and can be easily generalized to open universes pi[ . 

Since polarization fluctuations are typically at a part in a million, an order magnitude below the temperature 
fluctuations, to measure this signal requires high detector sensitivity, long integration time, and/or a large number 
of pixels. So far, only experimental upper limits have been obtained p^--|2^, with the current limit on the linear 
polarization being 16/iif |Q. Ground-based experiments being planned or built will probably achieve detection 
sensitivity using low-noise HEMT amplifiers as well as long hours of integration time per pixel. The MAP satellite 
will launch in 2000 and make polarization measurements of the whole sky in about 10^ pixels. If the polarization 
foreground can be successfully removed, MAP should marginally reach the detection level. For a detection of the 
magnetic polarization one would require either several years of MAP observations or the Planck mission . 
We expect that polarization measurements are as important as anisotropy in future missions. 

Previous full-sky studies of the polarization are mainly based on angular power spectrum estimators in Fourier 
space. Although electric- and magnetic-type scalar fields E and B in real space can be constructed, they must 
involve nonlocal derivatives of the Stokes components. In this paper, we will study in detail the auto- and cross- 
correlation functions of the Stokes parameters themselves in real space. Although the two approaches should be 
equivalent to each other, one can find individual advantages in different situations. We will follow the formalism of 
Ref. [|2|, expanding the Stokes parameters in terms of spin- weighted spherical harmonics. The expansion coefficients 
are rotationally invaria nt p ower spectra which will be evaluated using CMBFAST Boltzmann code developed by 
Seljak and Zaldarriaga ||2^. In Sec. El we briefly introduce the CMB Stokes parameters and their relation to spin- 
weighted spherical harmonics. Sec. ElD is devoted to discussions of the properties of the harmonics, the harmonics 
representation of rotation group, and the generalized addition theorem and recursion relation. In Sec. |^ we expand 
the Stokes parameters in spin-weighted harmonics, and briefly explain how to compute the power spectra induced by 
scalar and tensor perturbations. In Sec. ^ we derive window functions appropriate to detectors with Gaussian angular 
response in anisotropy and polarization experiments. The instrumental noi se of detectors in CMB measurements is 
treated in Sec. VI as white noise superposed upon the microwave sky. Sec. VII is to construct the auto- a nd cro ss- 
correlation function estimators corrected for noise bias in terms of the power spectra. As examples, in Sec. VIII we 
compute the means and variances of the estimators for different configurations of future space missions in standard 
cold dark matter models. Further, we outline the likelihood analysis of the experimental data in Sec. IX. Sec. ^ is 
our conclusions. 



II. STOKES PARAMETERS 



Polarized light is conventionally described in terms of the four Stokes parameters {I,Q,U,V), where / is the 
intensity, Q and U represent the linear polarization, and V describes the circular polarization. Each parameter is a 
function of the photon propogation direction fi. Let us define 

T^I-I (2.1) 

as the temperature fluctuation about the mean. Since circular polarization cannot be generated by Thomson scattering 
alone, V decouples from the other components. So, it suffices to consider only the Stokes components {T,Q,U) as 
far as CMB anisotropy and polarization is concerned. Traditionally, for radiation propagating radially along in 
the spherical coordinate system, see Fig. I, Q and U are deflned with respect to an orthonormal basis (d, b) on the 
sphere, which are related to {eg,e^) by 

d = 60, and b — —eg. (2.2) 

Then, Q is the difference in intensity polarized in the b and d directions, while U is the difference in the (d -I- 6)/\/2 
and (d — b)/\/2 directions p^. Under a left-handed rotation of the basis about through an angle ijj, 



a' \ _ / cos -0 — sin ip \ f a 
y ) I sinV' cosi/' ) \ b 



(2.3) 



or equivalently, 

i=(d' + .6') =e^'^i=(d + z6). (2.4) 
Under this transformation T and V are invariant while Q and U being transformed to [p7| 
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Q' \ _ f C0S2V' sin2V' \ f Q 
U' ) ^ [ -sin2V' C0S2V' ) \ U 



(2.5) 



which in complex form is 

Q'icr) ± iU'icr) = e^^'''' [Q{er) ± iU{er)] . (2.6) 

Hence, Q{n) ± iU{n) has spin-weii zM Therefore, we may expand each Stokes parameter in its appropriate 

spin- weighted spherical harmonics potp^ . 

Unfortunately, the convention in theory has a little difference from experimental practice. In CMB polarization 
measurements, usually the north celestial pole is chosen as the reference axis 63, and linear polarization at a point x 
on the celestial sphere is defined by 

Q{x) = Tn,s - Te,w, and U{x) = Tne,sw - Tnw,se, (2.7) 

where T^^s is the antenna temperature of radiation polarized along the north-south direction, and so on |2^ ]. In 
small-scale experiments covering only small patches of the sky, the geometry is essentially flat, so one can simply 
choose any local rectangular coordinates to define Q and U. Since an observation in direction x receives radiation 
with propagating direction h = — i, we have 

Q{x)^Q{h), and U{x) = -U{n). (2.8) 



III. SPIN- WEIGHTED SPHERICAL HARMONICS 

An explicit expression of spin-s spherical harmonics is f\ [^ , ^ 

'21 + 1 {i + my. {I -my: ' 



s ^ im ; ' 



(— l)™e*™'^ 



47r [i + sy. [i-sy 



sin2' 



where 

max(0, m — s) < r < min(Z — s, Z + to). (3-2) 
Note that the common spherical harmonics Yim = oi^m- They have the conjugation relation and parity relation: 

sYii,{e,cj,) = (-1)™+^ _,y,_„(0,0), (3.3) 

,y,™(7r-0,0 + 7r) = (-1)' ^sYi„,{e,^). (3.4) 
They satisfy the orthonormality condition and completeness relation: 

dn sYir„Ao, 0) sYUe, 4>) = SvAn'm. (3.5) 



sYii{e',(b') sYi,ni0A) ^ siq^' - cj))S{cose' - cose). (3.6) 



Im 

Therefore, a quantity 77 of spin-weight s defined on the sphere can be expanded in spin-s basis. 



'^Generally, a quantity -q will be said to have spin-weight s if it transforms as 77' = e^^^rj under the rotation ( [2.4| ) [p8[ . 
■^In Ref. psi, the sign ( — 1)™ is absent. We have added the sign in order to match the conventional definition for Yim- 
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r]{d,(j)) = ^riira sYln 



(3.7) 



where the expansion coefficients riim are scalars. 

The raising and lowering operators, d and d, acting on 77 of spin-weight s, are defined by |28 



drj = —(sin i 
dr] — — (sin^ 

When they act on the spin-s spherical harmonics, we have |p8| 



d a 

d .8 



{sin 9) "77, 
{siii9yr]. 



dsYlrn = [{I - S)il + S + 1)]^ s+lYl,n, 
'dsYhn = -[{l + - S + 1)]^ 

~l8sYi„, = -{l-s){l + s + l) sYi„,. 



(3.8) 
(3.9) 



(3.10) 

(3.11) 
(3.12) 



Using these raising and lowering operations, we obtain the generalized recursion relation for I — 2> max(|s|, |m|). 



/ + s 



{2l + l){2l - 1) 



cos6'5F/_i . 



{I + m){l ~ m) 
{21 + 1){1 + m - 1){1 - m - 1){1 - s - 1) 



+s 



{21 - 3){l + m){l - m){l + s - I) 
{2l + l){2l-l) 



sYl-2,r. 



{I + m){l - m){l - s){l + s ~ I) 



sin6' s_iy;_i, 



(3.13) 



This will be used for evaluating the correlation functions in Sec. |VII|. Table 1 lists explicit expressions for some low-^ 
spin-weighted harmonics, from which higher-/ ones can be constructed. 

The harmonics are related to the representation matrices of the 3-dimensional rotation group. If we define a rotation 
R{a, (3, 7) as being composed of a rotation a around 63, followed by /? around the new ej, and finally 7 around Cg , the 
rotation matrix of R will be given by psl 



An 



sYim{P,a)e' 



,21 + 1 

Let us consider a rotation group multiplication, 

R{a,(3,-j) = R{(t)',9\0)R-\^,9,0), 
where the angles are defined in Fig. |l|. In terms of rotation matrices, it becomes 

Di,s, (a, /3, -7) = E ^U(<^'' OMlmi^P, 0, 0), 



(3.14) 



(3.15) 



(3.16) 



which leads to the generalized addition theorem,^ 
E s,YC„,{9',cl,'),,Yi^{9,^)^ 



21 + 1 
47r 



{-ly^-^^ ^s,YuMc^)e^'''' 



(3.17) 



^This theorem was first derived in Eq. (7) of Ref. ||T^, which however does not give correct signs for the geometric phase 
angles, a and 7. Eq. (3.17) will be useful in the following sections. 
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IV. POWER SPECTRA 



Following the notations in Ref. ||l2|, we expand the Stokes parameters as 

T{h) = ^aT,imYi„i{n), 

hn 

Q{n) - iU{n) = '^a2jm 2Yi„i{n), 

hn 

Q{n) +iU{h) = ^a_2,im -2Yi„i{n). 



Im 



The conjugation relation (3.3) requires that 



2.lrn - (-l)"a2J- 



(4.1) 



(4.2) 



For Stokes parameters in CMB measurements, using the parity relation (3^), we have 

T{x) ^^{-iyaT,lrnYl„rix), 

Irn 

Q{x) + iU{x) = ^(-l)'a2j™ -2Yi^{x), 

hn 

Q{x) - iU{x) = ^(-l)'a_2.i™ 2Ylmix). 



(4.3) 



hn 



Isotropy in the mean guarantees the following ensemble averages: 

{o-T,l'm'(^T,lm) — CTlSl'l6m'm, 

{'^2Mm'(^2dm) — {C El + Cbi)Si'i6„ 

{(^2,l'm'(^~2,lm) — {C El — Cbi)Si'iS„ 

{°'T.l'm''^2,lm) = —CciSl'l6m'm- 



(4.4) 



Consider two points n'{0' , cj)') and n(6, (f)) on the sphere. Using the addition theorem ( ^.17 ) and Eq. (4.4), we obtain 
the correlation functions, 



o; I 1 

{T*{h')T{h)) = -i—CTiPi{cosP), 



{T*{n')[Q{n)+^U{n)])=-J2 



21 + 1 {1-2} 



47r 



(1 + 2) 



■CciPf {cos I3)e^ 



(4.5) 
(4.6) 
(4.7) 
(4.8) 



([Q(n') + iU{n')Y[Q{n) + iU(n)]) = ^ ^'^1±1{Cei + Cbi) 2li-2(/3, 0)e2'("-^), 
{[Q{n') - iU{n')nQ{n) + iU{n)]) = ^ ^'^I±1{Cei - Cbi) 2Yn{P, Q)e^'^^+'^\ 

where a, (i, and 7 are the angles defined in Fig. |^. Eq. ( |4.6| ) is the most general form of those found in Refs. . 
In the small-angle approximation, i.e. (3 « 1, a ~ 7, so Eq. (4.7) depe nds only on the separation angle (3. When n' 
and n lie on the same longitude, a = 7 = and hence E qs. (0 .4.7,4.8) depend only on (3. W hen n' and n lie on the 
same latitude, a + 7 = tt. Hence the phase angle in Eq. ( [4.8[ ) vanishes, and that in Eq. ( |4.7|) becomes equal to e"*'". 

A coordinate-independent set of correlation functions has been obtained by defining correlation functions of Stokes 
parameters {Qr,Ur) with respect to axes which are parallel and perpendicular to the great arc connecting the two 
points being correlated pit]. This prescription is indeed equivalent to the transformations: 



Qr{n') + iUr{h') = e-2n[Q(^') + iU{n')], 
Qr{n) + iUr{n) = e"2«"[Q(n) + iU{n)]. 



(4.9) 
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The authors in Ref. ||T^ expanded Qr and Ur in terms of tensor spherical harmonics. To calculate the two-point 
correlation functions between T, Qr, and Ur, they chose one point to be at the north pole and the other on the 
= longitude, and argued that the correlation functions depend only on the angular separation of the two points. 
Then, they had to evaluate the asymptotic forms for the tensor spherical harmonics at the north pole. Their results 
are simply equal t o the above correlation functions without the phase angles. Here, using the compact generalized 
addition theorem ( 3.17 ), we have given a general and straightforward way of obtaining the correlation functions. In 
addition, the phase information is retained. We will see in Sec. VI] and Sec. [X that these phase angles can be easily 
removed or evaluated in taking experimental data. 

Therefore, the statistics of the CMB anisotropy and polarization is fully described by four independent power spectra 
{Cti,Cei,Cbi,Cci) or their corresponding correlation functions. Here, we outline how to evaluate the spectra. The 
details can be found in Rcfs. f|14|Jlf: 



Since the four spectra are rotationally invariant, it suffices to consider the contribution from a single fc-mode of 
the perturbation, and then integrate over all the modes. In particular, the calculation will be greatly simplified if we 
choose k = 63. For scalar perturbations, the contribution of the fc-mode to (T{h), Q{n), U{n)) is (A^\ 0). For 
tensor perturbations, the contribution is 



(1 + 

V 



cos 



cos 20 A^^^ 
I cos 20 A^^^ 



(4.10) 



2 cos 6* sin 20 A^^ 



for -I—mode. The x-mode contribution is from making the replacements, cos 20 sin 20 and sin 20 — > — cos 20. The 
quantities A's are then computed by solving the Boltzmann hierarchy equations or by the line-of-sight integration 
method In the following sections, we will use CMB FAST Boltzmann code [|2^ to evaluate all Cxi's. 

V. WINDOW FUNCTION 

Due to the finite beam size of the antenna, any information on angular scales less than about the beam width is 
smeared out. This effect can be approximated by a Gaussian response function. 



aH{p, a) = — — e 



(5.1) 



where Ub, much less than 1, is the Gaussian beam width of the antenna, f3 and a are spherical polar angles with 
respect to a polar axis along the direction h{9,<j)). Therefore, a measurement can be represented as a convolution of 
the response function and the expected Stokes parameters. 



dR{l3,a)X{e' ,i 



(5.2) 



where X denotes T, or Q ± iU . This can be accounted by a mapping of the harmonics in Eq. (4.1) 



(5.3) 



From Eq. (3.17), we have 



47r 



21 



(5.4) 



Therefore, the convolution involves the integral, 



21 + 1 



(5.5) 



Making the approximation that a ~ 7 for df, << 1 and using the explicit expression (3.1), the integral ( |5.5| ) has a 
series solution as 
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(-lY{l-[(l-.s)il + s) + l]{^] + 



-a - s){i - s - + s){i + s - 1) 



{1(1 + I) -s')^ 



^—m.'s ■ 



~ (-l)^exp 
Hence, the mapping ( |5.3| ) is approximated by 

where is the window function, 

sWi = exp [- {1(1 + 1) - al] . 
When s = 0, it reduces to the usual window function in anisotropy case, 

oWi = Wi^e,^p[^l{l + l)ai]. 
The approximation ~ exp[— Z^ct^] works very weU for high Vs. 

VI. INSTRUMENTAL NOISE 



(5.6) 

(5.7) 
(5.8) 
(5.9) 



In the CMB experiment, a pixehzed map of the CMB smoothed with a Gaussian beam is created. In each pixel, 
the signal has a contribution from the CMB and from the instrumental noise. A convenient way of describing the 
amount of instrumental noise is to specify the rms noise per pixel tTpix, which depends on the detector sensitivity s 
and the time spent observing each pixel tpixi cTpix — s/y/t^. The noise in each pixel is uncorrelated with that in any 
other pixel, and is uncorrelated with the CMB component. Let Slpix be the solid angle subtended by a pixel. Usually, 
given a total observing time, tpix is directly proportional to flpi^. Thus, we can define a quantity w~^, the inverse 
statistical weights per unit solid angle, to measure the experimental sensitivity independent of pixel size ||3l|]: 



(6.1) 



Let us simulate the instrumental noise with a background of white noise superposed upon the microwave sky. The 
statistics of the white noise is completely determined by 

/ N * N \ _ — 1 r r 
\^T.,l'm'^T,l'm/ — Wj, OinOrn'm, 



N * N 



N * 

2.i'm'"-2,i' 



.) = o, 



(6.2) 



where the label N stands for noise, Wj,^ and Wp^ are constants to be dertermined. Then, the two-point correlation 
functions are 

21 + 1 



([Q^(n') + ^U''{h')nQ''{n) + ^U''{h)]) = ^ ^j^^^Wp' 2W1 ,Yi^,{P, 0)e'^^--^\ 



([Q^(n') - lU^inTiQ^'in) + ^U^in)]) = 



(6.3) 



Defining and be the rms anisotropy and polarization variances respectively, for small beam width we obtain 
from Eqs. (|U,U) that 
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. . . (6.4) 

Therefore, if we assume ilpix = 47rcr^, then the variances would be the pixel noise, and w^^p be the inverse statisical 
weights per unit soild angle: 

T.P T.P -1 r, f T,P\^ ,r. r-\ 

^ = %ix ' Wj^ p = f^pix (^CTpi^ j . (6.5) 

If both anisotropy and polarization are obtained from the same experiment by adding and subtracting the two 
orthogonal linear polarization states given equal integration times, then 

Kix)' = ^«ix)'- (6.6) 
If they are from different maps, the noise is uncorrelated. 



VII. FULL-SKY CORRELATION FUNCTION ESTIMATORS 



The CMB map is inevitably contaminated by instrumental noise and other known or unresolved foreground sources. 
However, the foreground contamination can be removed by observing the CMB at multi-frequencies and detecting its 
unique spectral dependence. After the removal of background contamination, the microwave map (denoted by label 
M) is made of the genuine CMB and instrumental noise: 



M 
^T,lm 



T,hm 



N 



(7.1) 



Thus the statistics of the noisy CMB map is induced from that of the CMB in Eq. [iA] and that of the noise in 



Eq. (6.2). Again note that the noise is uncorrelated with the CMB signal, i.e. (a^a) = 



Now we are going to construct the full-sky averaged correlation function estimators. Let us begin taking an average 
of a product of two spherical harmonics over the whole sky. 



{Y{;^,{n')YUn)}s = / dn' dnYi%^,{n')Yi^{n) 

Pi{cos (3)5i'i5 



in 



(7.2) 



where the curly bracke ts {} 5 denote a full-sky averaging at a fixed separation angle f3. The sky averaging can be 
done easily using Eqs. ( 3.17[ 3.5). We firstly transform Yi'm'{n') defined by a spherical coordinate system 63 to a new 
coordinate system n, and then performing azimuthal integration by rotating the transformed n' about n with a fixed 
separation angle (3. Finally, the remaining product of spherical harmonics with angle variables h is integrated over 
the whole sky. 

To generalize th e averaging procedure to spin-s spherical harmonics, some complications have to be taken. As we 
have seen in Eqs. (4^,4^,4^), multiplication of higher spin harmonics depend explicitly on local angles. Therefore, 
we define the full-sky averaging as 



1 



47r(2/ + 1) 



(7.3) 



Obviously, when si = S2 = 0, it reduces to Eq. (|7.2|). 

We define four full-sky averaged correlation function estimators. 
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Ct(/3) = {T^*(n')T^(n)}^ - ({r^*(n')r^(n)}^) 



Cc{0) = \ {r^*(n')[g^(n) +i[/^(n)] +h.c.}< 



= E \/^ - 2^p') 0), 
C_(/3) = i {[Q^(n') - iC/^(n')]*[g^(n) +i;7^(n)] +h.c.}^ 



(7.4) 



where 



nM _ 1 \ ' ^M* „M 

— 21 + 1 / ."'T,lm"'T,lmi 

m 

- ^ \^ n^* «^ 
~ 2Z + 1 ^ "'T,lm"'2,lm 



-CI 



pM 
'^±1 



21 + 



-E 
+ 1 ^ 



M* M 



(7.5) 



The ensemble mean of each estimator is 



(Ct(/3)) =^^CT,W^iP,(cos/3), 



4-77 



(7.6) 



And the covariance matrix can be constructed as 

^x'Y ^ {[CxiH') - {CximiCriP) - {Cvim) , (7.7) 
where X,Y = T,C,+,—. Here the prime denotes a different separation angle. The diagonal entries are given by 

^^'^ = ^ E(2« + 1) ('^T' + ^;(cos/3')P;(cos/3), 

Mc?'C = ^ E ^'^t^lllo!! + (^^' + "^t") {Cei + wp')] Wi 2M/,Pf(cos/3')Pf (cos/3), 



A-K^ An {l + 2)\ 
= ^ E ['^Ei + Cli + 2 {wp')' + 2 {Cei + Cm) Wp\ 



27r 
1 



2Wf 2li-2(/3',0) 2l;-2(/3,0), 



M_._ = ^ E [^Ei + Cli + 2 {wp')^ + 2 {Cei + Cbi) Wp'] -.W^ ^Yn{ii' , 0) ^Yn{(i, 0). 



(7.8) 



The ofF-diagonal entries are similarly calculated. In particular, the off-diagonal term of the submatrix M^'f {X, Y = 
+, -) is 



= ^ E [^Ei ~Cli + 2 {Cei - Cbi) wp'] 2Yi-i{0. 0) ^YniH, 0). (7.9) 
I 

In practical situations, a galaxy-cut on the CMB map is necessary due to radiation pollution along the galactic 
plane, and due to limited observation time, usually only a fraction of the sky would be sampled. For instance, the 
effective CMB coverage of the COBE DMR is 47r/, where / ~ 2/3. This incomplete sky coverage would generally 
induce a sample variance, whose size depends both on the experimental sampling strategy and the underlying power 
spectra of the fluctuations. It was found that the covariances calculated above scale roughly with sky coverage as /"^ 
for small-scale experiments . For large-scale experiments such as the COBE DMR, they scale roughly as 

0.446-I-0.542/-1 - 0.0079/"^ (7.10) 

valid for /^^ < 15 |^^. The difference from /^^ scaling is mainly due a large correlation angle in the large-scale 
experiment. 



VIII. CORRELATION MEASUREMENTS IN FUTURE MISSIONS 



The MAP and Planck missions plan to measure all-sky CMB anisotropy and polarization. It has been discussed how 
to construct the optimal estimators for the power spectra corrected for noise bias, and their corresponding variances 
from the all-sky map |l^,|l^,^. To estimate the level of signal and noise, we hereby give an alternative real-space 
analysis, evaluating the ensemble means and variances of the full-sky averaged correlation function estimators for the 
MAP and Planck configurations, i.e. 

Cx{e) = {Cx), ACx(e) = ML, (8.1) 

which are respectively given by Eq. ( |7.6| ), and Eq. ( |7.8| ) with 9' ~ 9. 

We assume the standard cold dark matter (sCDM) model: flo — 1, h ^ 0.5, VlBh? — 0.0125, and no reionization 
after the hydrogen recombination. Two extreme cases are evaluated: T/S = and T/S=l, where T and S are 
the anisotropy quadrupole moments induced respectively by tensor and scalar perturbations. All power spectra 
are computed by the CMBFAST code. The recursion relation ( 3.13D has been used for evaluating the spherical 



harmonics in the correlation functions. The results are plotted in Figs. ^[^, which are respectively Cx{9) attached 
with its variance ACx{9), where X = T,+, — , C. In making the plots, we have used beam width 9fwhm = 0.5°, 
where = 0.425 x 9fwhm- Typical values of the experimental sensitivity for MAP are w^^ = (0.1/iX)^ and 
Wp^ = (0.15/iif)^, while for Planck they are about a factor of 100 smaller. In Figs, ^j^, the thick and thin solid 
lines represent the cases with T/S ~ and T/S = 1 respectively. In each case, the ensemble average is denoted by a 
middle line sandwiched by two pairs of la lines. The outer pair of lines is for the MAP satellite while the inner pair 
for the Planck Surveyor. Note that in Fig. |^ the two pairs of la lines merge into a single pair, which means that the 
noise is dominated by cosmic variance rather than instrumental noise. 

The theoretical expectation of the rms polarization signal in sCDM models, [C+(0)]^/^, is at a level of IfiK. For the 
MAP experiment, the polarization signal to noise ratio S/N is about 1 — 2. The S/N ratio of the anisotropy-polarization 
correlation Cc{d) is about 3 — 4 at 6* ~ 1.3°, and the absence of tensor mode makes the cross-correlation significantly 
negative on few-degree scales. For Planck the corresponding S/N ratios are much higher. The MAP would likely 
detect the anisotropy-polarization correlation, which however is not sensitive to Cei or C'bi- The detection of the 
electric and magnetic components would require the Planck satellite. 

Another space misssion being planned is the Sky Polarization Observatory (SPOrt) on board the International Space 



Station during the early space station utilization period (2001-2004) |34|. The scope is to measure the polarization of 
the sky diffuse background radiation at an angular-scale of 7° for a large sky-coverage with four frequency channels 
between 20 GHz and 70 GHz. The experimental sensitivity is expected to be comparable to MAP. Again, we evaluated 
the ensemble means and variances of the full-sky averaged correlation functions, but in rcionized sCDM models with 
reionization redshifts Zri — 20 and 50. The results are plotted in Figs. ^-|^. The expected rms polarization S/N ^1 — 3 
for 20 < Zri < 50, while the anisotropy-polarization S/N ~ 1 — 2 at ~ 20°. 

A near-term, ground-based polarization experiment, called POLAR, is to measure CMB polarization at 7° scales 
for 36 pixels |35| ]. To reach a signal level of IfiK for a single pixel requires an integration time of about 120 hours and 
low-noise HEMT amplifiers of noise temperature of about lOK. The expected S/N ratio of the rms polarization is 
1 — 2 for reionized sCDM models with reionization redshifts 45 < Zri < 105 whereas the anisotropy-polarization 
correlation would be dominated by noise P-|Tl]] . 
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IX. LIKELIHOOD FUNCTIONS 



The most straightforward way to obtain the power spectra from the measured Stokes parameters is to perform a 
maximum Ukehhood analysis of the data. All of the information in the measurement is encoded in the likelihood 
function, which can properly take into account non-uniform detector noise, and sample variances. This is particularly 
an advantage for ground-based experiments which track tens or hundreds of spots in the sky to measure Q and 
U . The method offers a simple test of the consistency of the power spectra from map to map, and the correlation 
between maps. In fact, this method has been employed by the COBE/TfMR to determine the anisotropy quadrupole 
normalization from the two-point functions of the 4-year anisotropy maps containing about 4000 pixels [g^. However, 
as is known, for all-sky coverage in satellite experiments, especially small-scale measurements, the large amount of 
data involved in the computation makes the analysis inefficient. This problem may be reduced by using filtering and 
compression as in the case of anisotropy data. 

For a small number of measurements, such as the ongoing polarization experiment POLAR that measures Q and 
U by observing an annulus of regular 36 spots at constant declination [p5[ , the data set can be arranged as 



(9.1) 



where i,j — 1, .., 36. Since all data points lie on a same latitude, the expected theoretical correlation functions C±^j 
in this case are given respectively by Eqs. (4/7 4^ ) with a -I- 7 = tt, i.e. 



2/ + 1 
47r 



{CEl+CBl)2Wl2Yl-2{e^J,Q)e^ 



21 + 1 
47r 



(9.2) 



where 2W1 is the window function (5^) with a beam width appropriate to the experiment, 9ij is the separation angle 
between the ith and jth spots, and (which is a geometric function of Oij) is the angle between the longitude at 
the ith spot and the great arc connecting the ith and jth spots (refer to Fig. |l|). Thus we construct the likelihood 
function as 



C{Ce2,Cb2) — 



1 



VdetC 



exp 



-DC 

2 



(9.3) 



where the correlation matrix is 



(9.4) 



where Nij = 2{af)^6ij is the noise correlation matrix. The most likely electric and magnetic quadrupoles are then 
determined by maximizing the likelihood function over the theories. 

For all-sky measurements, the full likelihood function can be constructed as 



£{Ct2, Cc2, Ce2,Cb2) 



1 



VdetM 



exp 



(9.5) 



where M is the covariance matrix of the full-sky averaged correlation functions whose entries given by Eq. (7.7), and 
AC is a row vector with entries 



ACx{0) = Cx(^)mcasurcd — (Cx{d) ^-^^^^-^^i^ 



(9.6) 



where the first term is the two-point correlation function in the sky map obtained by performing the full-sky- 
averaging (7.3) of products of all map measured Stokes parameters with a fixed angular separation 6, and the second 
term is calculated from the ensemble mean of the corresponding operator listed in Eq. ( |7.4D without substracting 



off the noise, i.e. setting Wrp 



0. The tensor contribution can be analyzed by maximizing the likelihood 



fimction with the covariance submatrix Mx'y, where X,Y = The central value of Cb2 in a confidence- level 

plot significantly deviated from zero would indicate the presence of tensor mode. 
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X. CONCLUSIONS 



It is known that the two-point correlation functions of the Stokes parameters are exphcitly dependent on coordinates. 
A way of getting rid of this is to expand the Stokes parameters in terms of spin- weighted spherical harmonics, 
and to construct optimal angular power spectrum estimators. Although being useful for all-sky coverage satellite 
experiments, it is not suitable for near-term, ground-based polarization experiments. For a small number of observation 
points, the simplest way to compare data with theory is to perform a likelihood analysis with a correlation function 
matrix. Further, a likelihood analysis of a full-sky map using correlation functions is a challenge for development of 
computational algorithms. 

Several authors suggested to obtain coordinate-independent correlation functions by measuring Q and U with 
respect to axes which are parallel and perpendicular to the great arc connecting the two points being correlated . 
Here we gave the most general calculation of the two-point correlation functions of the Stokes parameters in terms 
of spin-weighted spherical harmonics, including the windowing function and instrumental noise. We obtained simple 
forms though they still explicitly depend on coordinates. However, the coordinate dependence can be eliminated 
by averaging over the whole sky, and the averaged correlation functions can be used to construct the covariance 
matrix in likelihood analysis of future CMB satellite data. Moreover, in ground-based polarization experiments, if a 
correct scanning topology is selected, the coordinate-dependence can be eliminated or simplified in some way and the 
correlation functions can be directly put in the correlation matrix of the likelihood function. 

Furthermore, we have calculated the signal to noise ratios from the two-point correlation functions for future 
anisotropy and polarization experiments. It is likely that MAP will detect the first anisotropy-polarization correlation 
signal. In fact, to complement the MAP measurement, a small- angle ground-based polarization experiment, targeting 
at a higher signal to noise for rms polarization, should be performed as to cross-correlate with the MAP high-precision 
anisotropy map. Surely, measurements of the microwave sky by Planck will push cosmology into a new epoch, as both 
CMB anisotropy and polarization can be precisely measured. On the other hand, the limit on Sunyaev-Zel'dovich 
distortion of Compton-y parameter, \y\ < 15 x 10~^, from FIRAS data [3^ constrains the reionization redshift 
Zri < 50 jSSj . This constraint is consistent with theoretical CDM model calculations , which predict an occurrence 
of reionization at a redshift 30 < Zj-i < 70, and most likely at Zj.; ~ 50. So, it is probable that SPOrt/ISS would 
observe a polarization signal. However, if large-scale polarization is not detected, then it would be an impact on 
cosmological theories. 
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ir2±2 = ±i T cos 0) Sin 6 e±2# 

i>3±2 = i^/lf [2sin3 6'-sin6l(l±cos6l)2] e±2# 

iYi±2 = ±^Y^[3sin6i(l ±cos6»)3 -5(lT5cos6')sin3 6'] e'^^''^ 

2y2±2 = iY^(lTcose)2e±2# 

2Y3±2 = ±j^y^[-(lTcos6l)3 + 5(lTcos6')sin2 6'] e^^i-/- 

2l4±2 = m\f^ [(It COS 6*)^ - 12(lTcos6i)2sin2 6i + 30sin^6l] e^^i-^^ 



TAB. 1. Some spin-weighted spherical harmonics with I = 2,3,4. 
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1 

FIG. 1. Spherical coordinates showing two unit vectors h'{9' ,4'') and h{9,4>) with separation angle (3. The angles between 
the great arc connecting the two points and the longitudes are 7 and a. 
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FIG. 2. Full-sky averaged correlation function and its variance, Grid) ± ACt(^), versus separation angle for MAP and 
Planck experiments. In each model, the correlation function is denoted by a middle line sandwiched by a pair of Icr lines, which 
represent the theoretical error expected in MAP and Planck experiments. This means that the error is dominated by cosmic 
variance rather than instrumental noise. 
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FIG. 3. Full-sky averaged correlation function C+{9) ± AC+{9). In each model, the correlation function is denoted by a 
middle line sandwiched by two pairs of la lines. The outer pair of lines is for the MAP configuration, while the inner pair for 
the Planck Surveyor. 



17 




18 




19 




-500 I ^ 
50 100 150 

9 [degrees] 



FIG. 6. Full-sky averaged correlation function and its variance, Ct{9) ± ACt(S), for SPOrt/ISS experiment. In each model, 
the correlation function is denoted by a middle line sandwiched by a pair of Icr lines. Solid and dashed lines denote respectively 
the models with reionization redshifts Zr\ = 50 and 20. 
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FIG. 9. Full-sky averaged correlation function Cc{0)±ACc{0) for SPOrt/ISS. 
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